Hawking radiation from a collapsing dust sphere and its back reaction at the 

event horizon — Weak value approach — 
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To see the back reaction of the Hawking radiation in a dynamical spacetime of the spherical 
gravitational collapse, we explicitly calculate the weak value of the energy-momentum tensor 
of the massless scalar field. The background geometry of a collapsing dust sphere is specified 



by using the Painleve-Gullstrand coordinates, in which the time coordinate coincides with 
the proper time of a free-falling observer and the metric tensor is regular at the event horizon. 
The result is that in the remote future the weak value diverges at the event horizon. We argue 
that since the semi-classical approximation of the Einstein equation in the sense of the weak 
value breaks down there, the future geometry of the spacetime cannot be the Schwarzschild 
geometry. 



PACS numbers: 03.65. Ta, 04.62,+v, 04.70.Dy 



I. INTRODUCTION 

In the context of classical general relativity, black holes only absorb particles but do not emit. 
, Absorption of particles or energy would make a black hole larger. Hawking mathematically showed 
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that the area A of the event horizon of each black hole does not decrease in time, 5A > 
classical gravity. Combining the black hole physics and thermodynamics, Beckenstein suggested 
that a black hole has its intrinsic entropy proportional to its surface area [2]. This property of 
black holes is analogous to the second law of thermodynamics. Bardeen, Carter, and Hawking have 
shown that classical black holes with the mass M, the angular momentum J and the electric charge 
Q satisfy the four laws in analogy to thermodynamics [3j|. The first law is that any two stationary 
axisymmetric black holes are related by 6M = (k/8tt)5A + H,n5J, where k is the surface gravity 
and is the angular velocity of the black hole. The surface gravity k of a stationary black hole is 
analogous to temperature in the zeroth law; k is constant over the event horizon. The statement 
of the third law is that it is impossible to reduce k to zero by a physical process {^J. 

However, the temperature of a black hole is to be absolute zero since if it were not zero a 
classical black hole would emit radiation. This puzzle was resolved in quantum field theory in 
curved spacetime by Hawking proving that a black hole formed by gravitational collapse will 



create and radiate particles at a steady rate, which is by now called the Hawking radiation [jj. 
The radiated particles are considered to have the thermal distribution with a finite temperature 
T = k/2tt, which is called the Hawking temperature. This result gives the explicit form of black 
hole entropy, S = A/ A. The area of the event horizon is likely to decrease due to the Hawking 
radiation. However, the sum of the quarter of the area A and the entropy of the radiation never 
decreases as the generalized second law claims [6]. Thus, the Hawking radiation conforms with 
black hole thermodynamics. 

Each different choice of the coordinate system would generally give a different physical phe- 
nomenon such as the Hawking radiation in quantum field theories, since different time slices give 
different vacuum states. The spacetime geometry of the spherically symmetric (uncharged) black 
hole and the empty region outside a spherical body is given by the Schwarzschild solution. The 
Schwarzschild metric can be represented not only in the Schwarzschild coordinates but also in 
the Eddington-Finkelstein ones, the Kruskal-Szekeres ones,_and so on Q|. Furthermore, it can 



The last one is particularly useful 



also be expressed in the Painleve-Gullstrand coordinates 
for exploring the black hole physics, since the Painleve-Gullstrand time coordinate coincides with 
the proper time of an observer freely falling from spatial infinity, while the Schwarzschild time 
coordinate coincides with the proper time of an observer being at rest at spatial infinity. The 
Painleve-Gullstrand metric tensor is spatially flat and has an off-diagonal element so that it is 
regular at the event horizon. The hypersurface of a constant time traverses the event horizon to 
reach the central singularity. Exploiting the convenient properties of the Painleve-Gullstrand coor- 
dinates, we constructed a simple expression for the gravitational collapse of a spherical dust star in 
the previous work [h]]. We introduced a generalized form of the Painleve-Gullstrand coordinates 
with the time coordinate being the proper time of a free-falling observer to describe both the inte- 
rior matter region and the exterior empty region of a collapsing sphere in a single coordinate patch. 
It describes the inside region of the event horizon as well as the outside. In a semi-classical theory 



using the 
e.g., in 



ainleve-Gullstrand coordinates, the calculation of the Hawking radiation is performed, 
111 ] as a tunneling process. 
In the present work, we are going to study the weak value of the energy-momentum tensor in 
the Painleve-Gullstrand coordinates. The weak value has been originally proposed by Aharonov 



and his collaborators 
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14i | in terms of the weak measurement. The weak value of an observable 
A is defined by < A w >:= (f\A\i) / {f\i), where \i) and (/| are an initial state and a final 
state, respectively. In general, the weak value is a complex number but can be measured by the 
weak measurement. Recently, many experiments have been performed to demonstrate the strange 
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values 
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1611. Since then the weak value has been recognized as a fundamental concept of quantum 



mechanics [171. Ilia]. Although it appeared in a different context, e.g., in quantum field theory, the 
explicit applications to black hole physics were pioneered by Englert and Spindel [lj]]. They 
replaced the conventional vacuum expectation value by the weak value of the energy-momentum 
tensor to examine the back reaction of the Hawking radiation in the semi-classical approximation. 
Although the weak value was not directly calculated, they argued that it is not regular at the future 
event horizon and therefore the semi-classical approximation may break down there. Namely, the 
Hawking radiation changes the structure of the classical event horizon completely. 

In this paper, to see the back reaction of the Hawking radiation we directly calculate the weak 
value of the energy-momentum tensor in a dynamical spacetime of the spherical gravitational 
collapse and check its regularity at the future event horizon. By a new method in which the two 
kinds of the Bogoliubov transformation are considered because of the distinction between solutions 
of the Klein-Gordon equation in the empty region and in the matter region, we demonstrate 
the Hawking radiation of scalar particles in the Painleve-Gullstrand coordinates. The Painleve- 
Gullstrand coordinate system that is simple and well-behaved at the event horizon makes a matrix 
element of the weak value directly calculable. As a result, we will see that in the remote future the 
weak value of the energy-momentum tensor certainly diverges at the event horizon. 

Traditionally people in quantum field theory have studied the vacuum expectation value of the 
energy-momentum tensor to study the back reaction of the Hawking radiation and found that the 
incoming energy flux is negative, which is consistent with the decrease of the black hole mass j^ol ]. 
We can reproduce this also in our weak value approach, since the vacuum expectation value is the 
weighted sum of the weak value with the weight being the probability to obtain the final state. Due 
to the large fluctuation it is obvious that the quantum state of matter cannot be semi-classical. 
However, we claim that the semi-classical gravity theory, in which the gravity field is treated 
as classical while the matter field is quantized and the energy-momentum tensor in the Einstein 
equation is replaced by its weak value, breaks down at the future event horizon. 

The organization of the present paper is as follows. We describe particle creation by the spherical 
dust collapse in the Painleve-Gullstrand coordinates in section [Hi To show the Hawking radiation 
and subsequent calculations, we specify the background geometry and the mode functions of the 
massless scalar field. The mode functions and the vacuum states reflect a physical feature of the 
coordinate system. In section IIIIl the weak value of the energy-momentum tensor of the field 
specified in section [H] is calculated as the back reaction of the Hawking radiation, which appears 
in the Einstein equation in the semi-classical approximation. Taking the advantage of the nice 
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properties of the Painleve-Gullstrand coordinates, we directly calculate the weak value at the 
event horizon to see its dependence on the radial coordinate. Section IIVI is devoted to summary 
and discussions. 



In this section we demonstrate the occurrence of particle creation in the case of the spherical 
dust collapse employing the Painleve-Gullstrand coordinates. The Painleve-Gullstrand coordinates 
in the Schwarzschild spacetime specify a distinct time slice from that of the Schwarzschild or the 
Kruskal-Szekeres coordinates. The Painleve-Gullstrand time coordinate has a physical meaning 
that it corresponds to the proper time of a free-falling observer not only in the empty Schwarzschild 
region but also in the collapsing dust region. In this coordinate system we show the particle creation 
in the black hole spacetime which emerges from the collapse. We ignore possible backscattering by 
the background geometry because it is not essential in the following discussion, and assume that if 
there exists particles which escape to the future null infinity they entirely come from the past null 
infinity propagating through the dust region (See FIG. [I]). 



II. PARTICLE CREATION BY A COLLAPSING MATTER 
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FIG. 1. The Penrose diagram for the emerging black hole. The in- vacuum state |in) is defined on the 
hypersurface of past infinity and the out- vacuum state |out) on the hypersurface of future infinity. 
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A. Background geometry and mode functions 



We assume that the geometry of the spherical collapsing dust spacetime is given by the metric 
form in the Painleve-Gullstrand coordinates, 

2 



-dt 2 + \ dr + 



2mj 



-dt + r 2 dfl 2 , 



(1) 



where mass functions are m+ = M (const.) in the exterior region of the dust sphere of mass M 
and m_ = 2r 3 /9t 2 in the interior region lfj|. In this case the dust sphere starts to collapse at rest 

1 /3 

in the remote past (t — > — oo) and the radius of its surface is given by r s (t) = (9M(—t) 2 /2j , 
(—oo < t < 0). The radius of the sphere is 2M at time t = — 4M/3 and is zero at t = 0. The 
metric above is of C 1 class in the whole spacetime. The relation between the Schwarzschild time 
coordinate t sc h and the Painleve-Gullstrand coordinates (t, r) is 

V2Mj 



dt S ch = dt 



1 - 2M/r 



dr, 



(2) 



and therefore, 



i sch = t - 2V2Mr + 2M log 



y/r/2M + l 



v / r/2M - 1 

The line element of the exterior geometry is rewritten as 



ds 



2M 

r 



1 dudv + r 2 dtt 2 , 



(3) 



(4) 



where r* = r + 2M\og\r/2M — 1| is the tortoise coordinate. The retarded and advanced null 
coordinates are respectively represented by 



u = i S ch - r* = t - £(r), 
u = tsch + r* = t + r?(r), 



where 



£(r) = r + 2\/2Mr + 4M log y / r/2M - 1 
^(r) = r - 2V2Mr + 4Mlog {y / r/2M + 1). 



(5) 
(6) 

(7) 
(8) 



In the asymptotic region r — )• oo, the solution of the massless Klein-Gordon equation □<& = is 
the linear combination of 



e- luJU Y lm ( 



and -e-T, m ( 



(9) 
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for the outgoing and incoming components, respectively. We define the in-modes (f>™ by the in- 
coming solution having the positive frequency with respect to the Painleve-Gullstrand time on the 
past null infinity , namely, 

1 



-e~ tuJV , u -> -oo, (10) 



where the spherical harmonics Yi m (6, <p) are suppressed (from now on we drop the angular part of 
modes). The modes </>™ are normalized with respect to the Klein-Gordon inner product 

(/«,/«') = -iJ d^V^g UMZ, - d^u*,) 



which now reads 



= -i Jdrr 2 f u lX-yJ^X)f^, (11) 



/+oo < > 

dv (rC) d v (rffi) • (12) 

The field operator is expanded in terms of the in-modes as 

d> = |^(a^C + «L nt C*)- (13) 

The in-vacuum state is defined as 

c£|in>=0 (14) 

by the annihilation operators a™ satisfying the commutation relation, [a™, a" 1 ,^] = 5{oj — a/). 

On the other hand, we define the out-modes ^>° ut as the outgoing solution on the future null 
infinity i.e., 

C^_^ e —«, w->+oc, (15) 
which is normalized with respect to the inner product 

/ + 00 < y 

(r4>T) ds (r^J U ) (16) 
-oo 

on the time slice of the remote future where its radial range is 2M < r < oo. The field can also be 
expanded as 

$ = Jdw (aT C Ut + C tf C ut *) (17) 
and the relevant vacuum state is the out-vacuum, 

a° ut |out) = 0. (18) 
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The particles propagating into the matter region from eventually reflect at the center of the 
dust sphere and go to the remote future. Focusing on the particles escaping through nearby the 
event horizon, it makes sense to take into account the intermediate modes on the interior matter 
region. The interior line element takes the form 



ds 2 _ = r 2 (-dT 2 + dR 2 ) + r 2 dn 2 = -r 2 dUdV + r 2 dft 2 , 



where the new time and radial coordinates are 

'6(-t)^ 1/3 



, R 



M J ' " \9M(-t) 2 
and then the null coordinates in the interior region are given by 



1/3 



(19) 



(20) 



U = T-R, V = T + R, 



(21) 



having the ranges — oo < U < —1 and — oo < V < 1. We denote the two independent exact 
solutions of the wave equation as 



out 



1 



-iuiU 



1 

7i 



1 + 



1 



e wV ( 1 + 



irujr 



iujT 
1 

iuiT 



(22) 
(23) 



which are the outgoing and incoming modes in the interior region, respectively. We may impose 
the boundary condition that the mode functions should be regular at the origin r = so that the 
intermediate functions in the matter region appear only in the combination ^ u = (ip™ 1 — ip 1 ^)- 

Massless particles propagate along null geodesies. Since the retarded and advanced null coor- 
dinates in the exterior geometry are continuously related to those defined in the interior at the 
surface of the dust sphere of radius r s , the two sets of coordinates have the functional relations 



u 



M 
~6~ 
M 
~6~ 



10 + 9C/ + C/ 3 -241og 
-10 + 91/ + T/ 3 + 241og 



-3- U 



2 

3 — V 



(24) 
(25) 



The retarded null relation (|24p diverges at the event horizon U = —3. Concerning the latest null 
ray that comes from and bounces off at the center of the sphere to escape to Eq. ([21 
reduces to 



U (u) ~ Uq — ae 



(26) 
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near the horizon u ~ +00 or U ~ Uq = —3, where k = 1/4M and a is a positive constant. On the 
other hand, Eq. (|25j) reduces for the latest incoming ray to 

V(v)~b + cv, (27) 

where 6 is a certain real number and c is a positive constant proportional to At. Therefore, the 
outgoing modes in the interior region may be regarded as the outgoing modes in the remote future 
(in the exterior region) 

VC* ~ —^ e - iu)U ^\ (28) 

and the incoming modes as the ones in the remote past 

!^ e -^W (29) 

In the interior region the intermediate mode solution appears only in the combination, ^ = 
(^° ut — if)™) by (I22p and (I23D . However, in the remote past there exists only the incoming modes 
ty w — > — tj)£ by (f29|) and in the remote future only the outgoing modes ty u — > ip° ut by ([28]) . 



B. Bogoliubov coefficients and Hawking radiation 

We calculate the Bogoliubov coefficients which relate the in- and out-modes. The out-modes 
can be expanded in terms of the in-modes as 

C ut = / doj" (cwC,, + /WC») (30) 

and also in terms of the intermediate modes as 

C lt = jdu'{A^u, + B^K>)- (31) 

Here we have neglected the backscattering by the exterior Schwarzschild geometry and assumed 
that all the particles reaching propagate through the interior region (See FIG.Q]). Particularly 
in the remote future Eq. (f3Tj) reduces to 

One can calculate the coefficients A UL0 > = (^SV'Jy*) an d Bom' = ~ (<#i ut > > usin g Ec l s - CSJ) 

and (|28p. to obtain 

B „„, j = " p [* (*? + u °)} r I 1 + "J • (33) 
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where T(z) is the gamma function. 

Comparing Eqs. (f30l) and (j3~Tj) . we see that the two kinds of the coefficients have a relation 



J dJ (4w ^ + B uu , %,) = J du" (cw C„ + /W <ffi) , (34) 
and on /" it reduces to 

dJ (-A^'C> - B^ifftr) = J du" (cwtffr/ + /WC») • (35) 

The coefficients a W£J // and /3 WW » which relate the in-modes </>™„ and the out-modes </>° ut can be 
obtained as 







do/{A^ (C'.C») + Aw (C*,C»)} ( 36 ) 



and 



/W'=/ ^'{^(C^C»)+^'(C'%C»)}- (37) 



o 



Calculating first the inner product in Eqs. (f36|) and (|37|l . one obtains 
and 

w,^- "^, ^-.-ft , /y. . (39) 

47rVa/w 77w + w "/c + ze 47r\/a;'w" ^ - W'/c - it 

Finally the integration over uj 1 in Eqs. (j36[) and (|37p gives the Bogoliubov coefficients which relate 
the in- and out-modes, 



aw 



, i (auj c) 1 
I = / — ex P 

R // I 27TVW(j/ / 



, vra; . „b-U 

± za; 

2k c 



r(l + ^Tia/'^°), (40) 



where T (z, t) is the incomplete gamma function. 

We express the out-modes in the remote future in the two different ways, i.e., by the in-modes 
in the remote past and by the intermediate modes in the matter region. The usual Bogoliubov 
coefficients are given by the integrals of the products of the Klein-Gordon inner products, coming 
from the relation of the two Bogoliubov transformations. Since the coefficients /3 WW » are nonzero, 
particle creation may occur in the situation considered in this section. 
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III. WEAK VALUE OF ENERGY-MOMENTUM TENSOR 

Imagine an observer in a small box which freely falls and crosses the horizon. He is supposed to 
quantum mechanically measure the Hawking radiation. Suppose that the result is zero particle or 
vacuum defined with respect to the outgoing modes. The initial state is the vacuum with respect 
to the incoming modes. What we want to know is the weak value of the energy-momentum tensor, 

,_ (out|T^|in) 
<J ^ >W '" (out|in) ' i4ij 

where |in) and (out| are the in- and out-vacuum states of quantum matter fields. 

In this section, we calculate the weak value of the energy-momentum tensor of the scalar field 
using the results in the previous section. The back reaction of the Hawking radiation is given as 
the weak value of the energy-momentum tensor of quantum matter fields. Note that the weak 
value naturally appears on the right-hand side of the Einstein equation, 

R»» - ^g^uR = 8vr < T Mi/ > w , (42) 

as justified, e.g., by the saddle point approximation in the path- integral quantization. 

In the remote future we have a reasonable form of the out-modes so that we can define the 
annihilation operators a° ut and therefore the out- vacuum state |out) by a° ut |out) = 0. Following 
the standard procedure of the Bogoliubov transformation, the in-vacuum state can be explicitly 
written in terms of the out-vacuum state as 



| in) = (out | in) exp 



|out) , (43) 



where V UUJ i = — (f3*a ,■ Then the weak value of the product of the fields is 

< > w = J du>dJ<fi?{x) (<W°' t V) + K^C'V)} > (44) 

where 5 WLJ i is a short hand of the Dirac delta function 5{ui — u/). The out-modes are given by the 
outgoing solution (|15p in the remote future. Their derivatives with respect to the time and radial 
coordinates are 

d t c ut = h^c*, (45) 

d^T = (i^ r e-r- 1 )C Ut - (46) 

Note that <9 r £ = (1 — \f2M/r)~ 1 = p~ l {r) in (j4Uj) diverges at the event horizon r = 2M. Using 
Eq. dHJ), the weak value of the energy-momentum tensor T^ v = d^<&d v <& — g flu d x $d\&/2 is given 
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by 

< T^(x) > w = J dudJ<^(x) {6^(l ) J t *(x)F^(x;u,u') - V^ J\x)G ixv {x;uj,u')} . (47) 

Each component of Fn V and contains the derivatives such as ()45|) and (|46p . Their behaviors 
near the event horizon are 

F rr , G rr = ujJp" 2 + 0(p- 1 ), (48) 

F tr , G tr = -ujJp- 1 + 0{ P °), (49) 

Ftt, Gu, and the other components = 0(p ). (50) 

The divergent behavior is a direct consequence of the radial derivative ()46|) of the outgoing solution, 
while < d x &d\Q > w is regular at the horizon. The leading behavior of the scalar < T^u^u" > w 
is p~ 2 , where = (u l ,u r ) = (1, —\j2Mjr) is the four-velocity of a radially free-falling observer. 

Therefore, we conclude that the weak value of the energy-momentum tensor diverges at the 
future event horizon. 



IV. SUMMARY AND DISCUSSIONS 



We have demonstrated the Hawking radiation of scalar particles in a dynamical spacetime of 
the spherical gravitational collapse and then calculated the weak value of the energy-momentum 
tensor to see the back reaction at the event horizon of the emerging black hole using the Painleve- 
Gullstrand coordinates in which the metric tensor of the background geometry is regular at the 
event horizon. It is shown that some components of the weak value of the energy-momentum 
tensor are not regular at the future event horizon due to the dependence of the mode functions on 
the radial coordinate. Therefore, the semi-classical approximation of the Einstein equation breaks 
down. 

In our paper we have mainly concerned with the weak value of the energy momentum tensor, 
while people normally studied the vacuum expectation value in the semi-classical approximation. 
We think the vacuum expectation value is not enough to see the behavior of quantum back reaction. 
For example, the semi-classical approximation breaks down, if the weak values for a particular final 
state is divergent there, while the vacuum expectation value, the average of the weak values, is 
finite there. Note the relation of the weak value < T^ u > w := (/| T^ u \i) / (f\i) with the expectation 
value (i\ T^ u \ i) : 

(i\T^\i) = J2\(f\i)\ 2 <T»v> w . (51) 
/ 
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Actually many calculations have shown that the vacuum expectation value is finite after regular- 
ization. We think that the divergences are cancelled out in the average. 

How can we physically interpret this? Consider an observer who is located far away from the 
star before its collapse. He will expect that a black hole will emerge and the Hawking radiation 
should occur. However, if someone detected a certain number of the Hawking particles or nothing 
at all, he should expect the resultant spacetime is far from the original Schwarzschild black hole 
because the conditioned value of the energy-momentum tensor is divergent at the future event 
horizon. The outcome spacetime is observer dependent. 

Englert and Spindel calculated the weak value in a clever but indirect way, e.g., via the trace 
anomaly. However, this obscures the physical origin of the divergence of the weak value. In our 
method, as mentioned above, it is obvious that the out-modes defined on the remote future give 
rise to the divergence at the event horizon. The out-mode function or rather the outgoing solution 
of the scalar field equation in the Schwarzschild spacetime breaks the semi-classical approximation. 
Hence we argue that the future geometry of the spacetime cannot be the Schwarzschild geometry. 

Several possible situations after a black hole evaporation have been considered by many people 
but with no consensus. Here we claim that the black hole formed by collapse never forms; the 
existence of the event horizon is inconsistent with the semi-classical Einstein equation. Even the 
lowest-order back reaction of the Hawking radiation breaks the structure of the event horizon 
critically. 

The Hawking radiation from black holes causes information loss paradox. Let us consider the 
situation in which a black bole is formed by a pure state. As the black hole emits thermal radiation 
and evaporates completely, the subsequent state is a mixed state. The transition from a pure state 
to a mixed state means a non-unitary evolution, which conflicts with quantum mechanics. More- 
over, the existence of the event horizon implies the existence of the domain of out-of-communication 
with outside of the horizon. The event horizon leads to the non-unitary evolution by itself. If one 
believes that the process of the evaporation of black holes is unitary, the event horizon seems to 
be never formed. 

We would like to briefly comment on possible effects on the back reaction of the Hawking 
radiation from the Planckian physics [21]. We have circumvented it in our modest claim that the 
Schwarzschild spacetime is not consistent with the semi-classical approximation. More precisely, 
the weak value of the energy-momentum tensor is described by the Painleve-Gullstrand coordinates 
which is regular at the future event horizon. Even if one tried the weak measurement there, this 
would not require energy of Planck scale. 
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